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Abstract  This  paper  investigates  the  effects  of  coflow  02  level  and  temperature  on  diffusion  flame  of  a  CH4/H2  jet 
in  hot  coflow  (JHC)  from  a  burner  system  similar  to  that  of  Dally  et  al.  The  coflow  02  mass  fraction  ( y*0  )  is  varied 
from  3%  to  80%  and  the  temperature  ( T*of )  from  1200  K  to  1700  K.  The  Eddy  Dissipation  Concept  (EDC)  model 
with  detailed  reaction  mechanisms  GRI-Mech  3.0  is  used  for  all  simulations.  To  validate  the  modeling,  several  JHC 
flames  are  predicted  under  the  experimental  conditions  of  Dally  et  al.  [Proc.  Combust.  Inst.,  29  (1),  1147-1154 
(2002)]  and  the  results  obtained  match  well  with  the  measurements.  Results  demonstrate  that,  when  y*0  decreased, 
the  diffusion  combustion  is  likely  to  transform  from  traditional  combustion  to  MILD  (Moderate  or  Intense 
Low-oxygen  Dilution)  combustion  mode.  When  T*of  is  higher,  the  temperature  distribution  over  the  whole  domain 
trends  to  be  more  uniform.  Reducing  or  Tc*of  leads  to  less  production  of  intermediate  species  OH  and  CO.  It  is 
worth  noting  that  if  yf  is  high  enough  ( y*0  >80%),  increasing  y*0  does  not  cause  obvious  temperature  increase. 
Keywords  jet  in  hot  coflow,  moderate  and  intense  low-oxygen  dilution  combustion,  diffusion  flame,  intermediate  specie 


1  INTRODUCTION 


Moderate  or  Intense  Low-oxygen  Dilution  (MILD) 
combustion  belongs  to  the  volume  combustion  which 
is  full  of  the  furnace.  Its  general  characteristics  are 
low  reaction  rate,  uniform  thermal  distribution,  low 
combustion  temperature,  low  noise  and  extremely  low 
combustion-generated  pollutants  [1-5].  Compared  to 
traditional  combustion,  it  can  improve  thermal  effi¬ 
ciency  (thus  reducing  CO2  emission)  by  more  than 
30%  and  reducing  NOx  emission  by  more  than  70%  [5]. 
This  combustion,  debuting  in  early  1990s,  has  been 
successfully  utilized,  particularly,  in  steel  and  metal¬ 
lurgy  industries  of  many  countries,  e.g .,  Germany, 
Italy,  Japan,  United  States,  Sweden  and  China.  Due  to 
its  high  thermal  efficiency  and  low  pollution  emission, 
MILD  combustion  is  considered  as  one  of  the  most 
promising  combustion  technologies  in  the  2 1  st  century. 

In  the  late  two  decades  or  so,  MILD  combustion 
has  been  increasingly  studied.  Szego  et  al.  [6]  studied 
the  effects  of  the  inlet  conditions  (i.e.  pre-heating 
temperature  of  the  air,  equivalence  ratio,  heat  dissipa¬ 
tion  rate  of  the  furnace,  fuel  dilution  level)  on  estab¬ 
lishing  MILD  combustion  on  a  laboratory-scale  fur¬ 
nace.  They  investigated  the  CO,  NOx  emission  and 
stability  characteristics  of  a  parallel  jet  combustion 
system.  Mi  et  al.  [7,  8]  performed  experiments  and 
computation  on  the  same  furnace  to  study  the  effects 
of  the  inlet  momentum  of  the  premixed  jet  of  fuel  and 
air  on  MILD  combustion  and  emissions  (i.e.  CO,  NOx). 
They  suggested  that  there  was  a  critical  momentum 


rate  of  fuel-air  mixture  above  which  MILD  combustion 
could  be  realized.  Li  et  al.  [9]  studied  the  effect  of 
pre-heating  temperature  of  air  on  the  characteristics  of 
CH4  combustion  and  NO  emission  by  computation. 
They  found  that  the  NO  emission  concentration  in¬ 
creased  exponentially  with  pre-heating  air  temperature. 
Li  et  al.  [10-12]  and  Mi  et  al.  [13,  14]  performed  ex¬ 
periments  and  computations  to  study  effects  of  inlet 
conditions  of  premixed  and  non-premixed  jets  on  re¬ 
alizing  MILD  combustion.  They  reported  that  when 
the  equivalence  ratio  was  near  1,  the  combustion 
mode  was  still  in  the  MILD  mode  and  NOx  emissions 
were  low  but  extremely  high  CO  and  H2  emissions 
were  obtained.  The  studies  [10-14]  focused  on  the  ef¬ 
fects  of  inlet  conditions  but  not  on  those  of  local 
quantities  of  the  flue  gases  inside  furnaces.  This  is  not 
sufficient  for  the  understanding  of  the  mechanism  of 
MILD  combustion. 

In  the  JHC  open  flame  device  used  by  Dally  et  al. 
[15],  temperature  and  02  level  of  the  surrounding  hot 
coflow  can  be  well  controlled.  Thus,  the  effects  of 
thermal  field  and  oxygen  dilution  of  the  recirculated 
hot  flue  gases  in  furnaces  may  be  simulated  by  using 
the  JHC  device  after  removing  the  cold  air  coflow. 
Many  simulations  by  computational  fluid  dynamics 
(CFD)  of  the  MILD  combustion  [7-11,  16-19]  have 
been  performed.  Christo  et  al.  [16]  and  Mardani  et  al. 
[17]  simulated  the  JHC  experiment  of  Dally  et  al  [15]. 
Comparison  between  the  computation  and  the  experi¬ 
ment  showed  that  the  Eddy  Dissipation  Concept  (EDC) 
model  [19]  and  detailed  mechanisms  by  GRI-Mech  3.0 
[20]  could  well  predict  the  characteristics  of  flow  and 
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flame  in  the  MILD  combustion.  These  well-proved 
models  are  used  in  the  present  study. 

The  present  work  is  to  first  verify  the  RANS 
modeling  by  simulating  the  JHC  combustion  under  the 
same  condition  of  Dally  et  al  [15].  Then,  an  investiga¬ 
tion  is  carried  out  for  the  effects  of  the  coflow  O2  level 
( y*0  )  and  temperature  ( T*0 f )  on  the  JHC  diffusion 
flame  after  removing  the  outer  cold  air.  The  objective 
is  to  clarify  how  y*Q  and  7c*f  of  the  hot  flue  gas  in 
furnaces  impact  the  MILD  combustion  of  gaseous  fuel. 

2  COMPUTATIONAL  DETAILS 

2.1  Configuration  of  the  JHC  burner  of  simulation 

Figure  1  shows  the  structure  of  the  JHC  burner 
[15],  from  which  the  modeled  flames  are  produced.  It 
consists  of  an  insulated  and  cooled  central  jet  (inner 
diameter  D  =  4.25  mm)  and  an  annulus  nozzle  (inner 
diameter  =  82  mm)  with  an  internal  burner  mounted 
upstream  of  the  exit  plane.  The  internal  burner  provides 
hot  combustion  products  which  are  mixed  with  air  and 
nitrogen  via  oxidant  inlets  to  control  the  temperature 


and  O2  levels  in  the  mixture.  The  fuel  is  a  mixture  of 
H2  and  CH4,  equal  in  volume.  In  the  experiment  of 
Dally  et  al.  [15],  the  JHC  burner  was  mounted  on  a 
wind  tunnel  which  provided  a  co-axial  surrounding  air 
stream  at  room  temperature  and  with  the  same  veloc¬ 
ity  as  the  hot  coflow  (3.2  nvs-1).  Their  mass  fraction 
of  O2  in  the  hot  coflow  (yo0 )  was  operated  at  3%,  6% 
and  9%. 

To  validate  the  RANS  modeling,  the  flames  from 
the  above  JHC  burner  are  simulated  under  the  identi¬ 
cal  conditions  of  Dally  et  al.  [15],  as  Table  1  shows. 
However,  to  study  the  effects  of  coflow  oxygen  and 
temperature,  the  boundary  conditions  of  the  JHC 
burner  [15]  is  modified  by  replacing  the  outer  air 
stream  with  all  hot  coflow,  as  Fig.  2  (b)  shows.  The 
removal  of  the  cold  air  is  to  ensure  that  y*Q  and  7c*f 
remain  the  same  in  the  whole  computational  domain. 
This  modification  makes  it  possible  to  better  investi¬ 
gate  the  effects  of  the  coflow  y*0  ,  7c*f  in  the  whole 
computational  domain  and  mimic  the  effects  of  the 
local  O2  level  and  temperature  of  the  hot  recirculation 
flue  gas  inside  furnaces  with  no  interference  by  the 
outer  cold  air. 


computational  domain 


oxidant 


perforated  plate 

X 


cooling  gas  inlet 

cooling  gas  outlet 


!  internal  burner 


A  1 

1 

oxidant 


520 


>  hot  coflow 

> 

fuel 


^  hot  coflow 


20 

:h« — >\ 


<N 

00 


'I 

rl  - 
I' 

/ 


o 

o 

o 


Figure  1  Structure  of  the  JHC  burner  (in  unit  of  mm) 


Table  1  Inlet  conditions  of  Dally  et  a/.’s  experiment  [15]  for  firing  CH4/H2 
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(a)  Experimental  cases  of  Table  1  [12]  (in  unit  of  mm) 
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(b)  Modified  JHC  burner  for  CFD  cases  of  Table  2  (in  unit  of  mm) 


(c)  The  structure  grid 


Figure  2  Boundary  conditions  (a)  experimental  cases  of  Table  1  [15]  for  validation  of  numerical  program,  (b)  modified  JHC 
burner  for  CFD  parametric  study  in  Table  2,  and  (c)  the  structured  grid 


2.2  Computational  models 

Due  to  the  symmetry  of  the  system,  a  geometri¬ 
cally  simplified  axisymmetric  computational  model 
(2D)  is  constructed,  as  shown  in  Fig.  2.  A  primary 
orthogonal  structured  mesh  with  about  30000  cells  in 
Fig.  2  (c)  is  used  to  simulate  flames  after  verifying  the 
grid-independency  of  the  results  using  a  finer  grid  with 
130000  cells.  The  computational  results  for  different 
grids  are  plotted  by  the  dimensionless  radial  location 
r/D  on  the  horizontal  axis,  the  temperature  (T)  on  the 
left  vertical  axis  and  OH  mass  fraction  (you)  on  the 
right  vertical  axis,  see  Fig.  3.  Results  (i.e.  temperature, 
OH  mass  fractions)  keep  in  high  consistency  for  dif¬ 
ferent  grids,  as  Fig.  3  shows.  The  origin  of  the  coordi¬ 
nate  system  is  at  the  fuel  nozzle  exit  (Fig.  2).  The  right 


and  upper  boundary  conditions  are  set  to  the  pressure 
outlets,  namely,  the  specification  of  a  static  (gauge) 
pressure  enforced  at  the  outlet  boundary.  The  value  of 
the  specified  static  pressure  is  used  only  when  the 
flow  is  subsonic.  All  other  flow  quantities  are  extrapo¬ 
lated  from  the  flow  in  the  interior. 

In  the  present  simulations,  the  standard  k-s  model 
with  the  standard  wall  functions  is  used  as  the  turbulent 
model.  The  standard  k-s  model  [21]  is  a  semi-empirical 
model  based  on  model  transport  equations  for  the  tur¬ 
bulence  kinetic  energy  (k)  and  its  dissipation  rate  (s): 
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Figure  3  Temperature  and  OH  mass  fraction  distribution  at  x  =  30  mm  for  the  two  different  grids  (curves  for  the  gird  with 
about  30000  cells,  and  symbols  for  the  gird  with  about  130000  cells) 


dips)  d^psiij)  d  \(  pt  \  ds 

— — — +  — - — l—  = -  //  +  — - + 

dt  dxf  dxj  ^  cik  J  dxj 

Cl£j{Gk+C,£Gb)-C2£p^-  +  S£  (2) 

In  these  equations,  p  is  the  fluid  density,  xt  and  Xj  (/, 
j  -  1,  2,  3)  are  the  coordinates,  ut  is  the  velocity  com¬ 
ponents  in  three  directions,  p  is  the  dynamic  viscosity, 
pt  is  the  turbulent  viscosity.  G/v  represents  the  genera¬ 
tion  of  turbulence  kinetic  energy  due  to  the  mean  ve¬ 
locity  gradients  whereas  Gb  is  the  generation  of  turbu¬ 
lence  kinetic  energy  due  to  buoyancy.  YM  represents 
the  contribution  of  the  fluctuating  dilatation  in  com¬ 
pressible  turbulence  to  the  overall  dissipation  rate.  C\£ , 
C2s  and  Cy-  are  constants,  oy  and  o£  are  the  turbulent 
Prandtl  numbers  for  k  and  £,  respectively.  Sk  and  Se  are 
user-defined  source  terms.  The  pt  is  computed  by 
combining  k  and  s  as  follows: 

M,  =  PC M  4-  (3) 

k 

where  C,(  is  a  constant.  The  model  constants  C\e,  C2s, 
Cfh  Ok  and  o£  are  usually  given  as  follows  [21]: 
C\£—  1.44,  C2s=  1.92,  C/f=0.09,  Ok=  1.0  and  o£  -  1.3. 
As  suggested  in  Ref.  [16],  the  coefficient  C\£  in  the 
eddy  dissipation  equation  should  be  varied  from  the 
standard  value  of  1.44  to  the  value  of  1.6,  which  has 
been  found  to  obtain  the  best  agreements  with  ex¬ 
periments.  The  present  study  has  therefore  taken  that 
C\£=  1.6.  The  discrete  ordinate  (DO)  radiation  model 
[18]  with  weighted  sum  of  gray  gas  model  (WSGGM) 
is  employed  as  the  radiation  model. 

The  EDC  model  [19]  with  the  detailed  chemical 
kinetic  mechanisms  (GRI-Mech  3.0  [20])  is  applied  for 
the  modeling  of  reactions.  In  the  present  study,  the 
NOx  emission  characteristics  are  not  investigated. 
Thus  the  species  and  reactions  involved  with  element 
N  included  in  GRI-Mech  3.0  are  eliminated.  A  modi¬ 
fied  GRI-Mech  3.0  which  consists  of  36  species  and  a 
total  of  219  reversible  reactions  is  applied.  This  has 
also  been  proved  in  Refs.  [22,  23].  In  EDC  model,  the 


species  conservation  equation  for  chemical  species 
takes  the  following  general  form: 

^l  +  VipvYi)  =  -VJi+Ri  (4) 

ct 

where  7/  is  the  local  mass  fraction  of  each  specie  /,  v  is 
the  velocity  vector,  Jt  is  the  diffusion  flux  of  specie  /, 
and  Rj  is  the  net  rate  of  production  of  species  i  by 
chemical  reaction.  The  fine  length  scale  (c)  and  the 
residence  chemical  time  scale  (r)  of  fluid  in  the  fine 
flow  structures  is  expressed  by 

(ve\  fvY/2 

#  =  Q  —  ,  r  =  CT  -  (5) 

\k  J  \£j 

where  and  CT  are  time  scale  constants  equal  to 
2.138  and  0.408,  respectively.  Rt  is  computed  for  each 
mean  species  as 

Rt=  (Y?-Yt)  (6) 

t( i-^y 

where  Y*  is  the  fine-scale  species  mass  fraction  after 
reacting  over  r.  The  evolution  of  Y*  depends  also  on 
the  chemical  kinetic  mechanism.  To  reduce  the  com¬ 
putational  cost  of  time  integration,  the  in-situ  adaptive 
tabulation  (ISAT)  model  of  Pope  [24]  is  used.  The  ac¬ 
curacy  is  checked  by  lowering  the  ISAT  error  toler¬ 
ance  and  ensuring  results  having  no  further  change. 
The  specific  heat  of  species  (Cp)  is  set  to  change  with 
the  temperature. 

The  SIMPLE  algorithm  method  is  utilized  to 
solve  the  pressure-velocity  coupling.  The  second-order 
upwind  scheme  is  employed  for  discretizing  the  equa¬ 
tions  in  order  to  improve  the  accuracy  of  the  simula¬ 
tions.  Convergence  is  obtained  when  (1)  the  residuals 
are  less  than  10~6  for  the  energy  and  10~5  for  all  the 
other  variables  and  (2)  the  variations  of  the  down¬ 
stream  outlet  temperature  and  velocity  are  allowed  to 
be  within  1  K  (<0.1%  of  the  outlet  average  temperature) 
and  0.001  nvs-1  (<0.1%  of  the  outlet  average  velocity), 
respectively. 
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Figure  4  Comparison  of  numerical  and  exprimental  profiles  [15]  of  the  temperatures  (7)  and  other  species  mass  fractions 
(  Joh  >  y0l  ’  yco  )  at  *  =  30  mm 
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3  RESULTS  AND  DISCUSSION 

3.1  Comparison  between  the  present  computation 
and  previous  experiment  [15] 

To  verify  the  numerical  predictions,  comparison 
is  made  between  the  numerical  results  and  the  ex¬ 
perimental  data  [15].  Fig.  4  compares  the  numerical 
and  exprimental  profiles  of  the  temperatures  ( T)  and 
species  mass  fractions  ( y0H  ,  y0 ^  ,  yco )  obtained  at 
x=  30  mm.  It  is  shown  that  the  radial  profiles  of  the 
temperature  and  OH  mass  fraction  are  predicted  quite 
satisfactorily.  The  difference  of  the  temperature  peak 
between  the  computation  and  the  experiment  is  smaller 
than  1.5%  (about  20  K).  Fig.  4  (b)  demonstrates  that 
the  variation  trends  of  the  CO  and  O2  mass  fractions 
agree  reasonably  well  with  the  experimental  data. 
Fig.  5  compares  the  present  predictions  (curves)  and 
the  previous  measurements  (dots)  [15]  of  the  mixture 
fraction  (£)  obtained  from  Bilger’s  formula  [25]  for 
y*0  =3%,  along  the  jet  axis  and  the  radial  direction  at 
x  =  30  mm,  60  mm  and  120  mm.  Obviously,  the 
predictions  agree  quite  well  with  the  measurements. 
This  suggests  that  the  RANS  modeling  used  presently 
can  approximately  capture  the  features  of  the  diffusion 
combustion  process. 

3.2  Effects  of  outer  cold  air  on  the  flame  structure 

Figure  6  shows  the  flames  of  Dally  et  al.  [15]  for 
y*0  =  3%,  6%  and  9%.  Evidently,  when  y q  =3%, 
there  is  no  luminous  flame  in  the  region  near  the  exit 
of  the  fuel  jet;  namely,  it  is  in  the  MILD  combustion 
mode  over  there.  However,  at  about  100  mm  above 
the  jet  exit,  the  outer  cold  air  is  entrained  into  the  hot 
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Figure  5  Comparison  of  the  present  predictions  (curves) 
and  the  previous  experiments  (dots)  [15]  of  the  mixture 
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coflow  and  starts  to  influence  the  flame,  consequently 
making  it  luminous  due  to  higher  oxidation  rates.  In 
the  experiment,  only  the  data  in  the  region  near  the 
exit  of  fuel  jet  (i.e.  x^30  mm)  may  be  appropriately 
used  to  investigate  effects  of  y*0  or  Tc*  f  on  the  MILD 
diffusion  combustion. 

Temperature  contours  of  the  JHC  flame  for  the 
case  of  Dally  et  al.  [15]  and  that  without  an  outer  cold 
air  coflow  are  shown  in  Figs.  7  (a)  and  7  (b),  respec¬ 
tively.  As  expected,  the  temperature  distributions  for 
the  two  cases  are  very  different  in  the  downstream 
region  at  x>30  mm.  In  the  case  with  an  outer  cold  air 
coflow,  similar  to  the  experiment,  the  temperature  in  the 
mixing  layer  seems  to  be  much  higher  than  7C*  f  down¬ 
stream  fromx=  100  mm,  due  to  more  oxygen  entrained 
into  the  hot  coflow.  By  comparison,  see  Fig.  7  (b),  the 
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Figure  6  Photographs  of  flames  reproduced  from  Dally  et 
al .  [15] 


temperature  is  only  slightly  higher  than  7c*f  in  the 
mixing  layer  of  the  jet  with  no  cold  air  coflow.  Obvi¬ 
ously  the  outer  cold  air  has  a  significant  impact  on  the 
downstream  combustion  process.  Thus,  the  following 
study  is  focused  on  the  cases  without  air  outside,  to 
investigate  the  effects  of  y q2  and  T  cof. 

3.3  Effects  of  coflow  02  level  (  Jq  )  on  the  flame 
structure 

As  listed  in  Table  2,  y*0  is  varied  from  80%  to 
3%  with  7c*f  kept  at  1300  K.  Fig.  8  shows  that  as  jPq 
decreases,  the  peak  temperature,  OH  and  H2CO  mass 
fractions  decrease.  The  temperature  and  OH  fields 
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Temperature  contours  of  the  JHC  flame  at  y(*}  =  3%  [(a)  with  an  outer  cold  air  coflow;  (b)  without  an  outer  cold  air  coflow] 


Table  2  Specifications  of  numerical  conditions  for  the  present  research 


Fuel  jet  (CH4/H2) 

Oxidant  coflow  (mass  fraction) 

Case 

Re 

T*/  K 

C  /K 

* 

v% 

*  /o/ 

y-nj/o 

*  10/ 
Fco2//o 

varying  y*2 

1 

9482 

305 

1300 

80 

8 

6.5 

5.5 

2 

9482 

305 

1300 

60 

28 

6.5 

5.5 

3 

9482 

305 

1300 

40 

48 

6.5 

5.5 

4 

9482 

305 

1300 

23.3 

64.7 

6.5 

5.5 

5 

9482 

305 

1300 

12 

76 

6.5 

5.5 

6 

9482 

305 

1300 

6 

82 

6.5 

5.5 

7 

9482 

305 

1300 

3 

85 

6.5 

5.5 

varying  T’o{ 

8 

9482 

305 

1700 

6 

82 

6.5 

5.5 

9 

9482 

305 

1600 

6 

82 

6.5 

5.5 

10 

9482 

305 

1500 

6 

82 

6.5 

5.5 

11 

9482 

305 

1400 

6 

82 

6.5 

5.5 

12 

9482 

305 

1300 

6 

82 

6.5 

5.5 

13 

9482 

305 

1200 

6 

82 

6.5 

5.5 
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Figure  8  Temperature  field,  OH  distribution  and  H2CO  distribution  in  the  whole  computational  domain  with  different 


Figure  9 
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CH4  distribution  and  CO  distribution  in  the  whole  computational  domain  with  different  Jq 


become  more  uniform.  OH  exists  primarily  in  the 
high-temperature  area  while  H2CO  is  mainly  in  the 
low-temperature  region.  H2CO  is  an  intermediate  species 
in  the  oxidation  of  most  hydrocarbon  fuels.  It  is  found 


typically  in  MILD  combustion  regime  [26].  There  is  an 
area  near  the  fuel  jet  exit  where  the  H2CO  mass  fraction 
is  0.  It  means  no  reaction  occurs  in  that  area.  Fig.  9 
shows  that,  as  y q2  decreases,  CH4  and  CO  mass 
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Figure  10  Radial  profiles  of  mean  temperature  for  Cases 
1-7 
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Figure  11  Radial  profiles  of  mean  mass  fraction  of  OH 
for  Cases  1-7 
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fraction  decrease.  That  is,  decreasing  4  generally 
leads  to  less  consumption  of  CH4,  thus  less  production 
of  CO.  In  addition,  it  takes  a  longer  distance  to  bum 
out  CH4  and  CO  when  4  decreases. 

Additionally,  Dally  et  al.  [15]  showed  that,  in  the 
case  with  4  =  3%,  at  x  =  30  mm,  the  temperature  rise 
is  about  100  K,  OH  mass  fraction  peak  is  about 
3.5 xlO"4  and  no  luminous  flame  can  be  seen  in  their 
photographs.  Namely,  the  combustion  is  in  the  MILD 


mode.  When  y q  increases  to  9%,  they  found  that  at 
x  =  30  mm,  the  peak  temperature  is  about  400  K  higher 
than  the  7c*f  and  OH  mass  fraction  peak  is  about 
13.2><10~4,  thus  traditional  luminous  flame  occurs. 
Thus,  in  the  present  study,  the  experimental  data  is 
used  to  judge  the  MILD  combustion  obtained  when  (1) 
the  temperature  rise  is  smaller  than  300  K  and  (2)  OH 
mass  fraction  peak  is  smaller  than  6xl0"4.  As  Figs. 
10-12  show,  it  can  be  easily  inferred  that  Case  7  is  in 
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the  MILD  mode  while  Cases  1~6  in  traditional  mode. 
Moreover,  in  the  MILD  mode  in  Case  7,  the  maxi- 
mums  of  CH3,  H2CO  and  CO  (the  unbumt  species) 
mass  fraction  are  1.1  X1(T4,  1.5 X1(T4  and  4.13xl0-3, 
respectively.  As  y*0  decreases,  classical  flame  mode 
tends  to  be  transformed  to  MILD  one. 

As  demonstrated  in  Fig.  12,  the  temperature  rise 
( A Tmax  =  rmax  -  T*0 f )  over  the  whole  computational 
domain  increases  with  y q  and  when  y q  is  high 
enough  ( y*0  >80%),  the  temperature  rises  to  a  con¬ 
stant.  Namely,  increasing  y q  cannot  improve  the 
temperature  rise  infinitely. 


for  Cases  1-7.  OH  is  a  vital  intermediate  and  ex¬ 
tremely  useful  in  diagnostics  usually  as  flame  marker 

[27].  As  demonstrated,  decreasing  y q  causes  rT 
and  r( 

a  little  larger  than  rT  (about  2  mm).  This  could  be 

1  max 

well  explained.  The  magnitude  of  OH  mass  fraction  is 
mainly  governed  by  two  elementary  reactions  as  fol¬ 
lows  [28]: 


0H  increase.  Worth  noting  that  r0H  is  always 


h+o2  — 

— »OH  +  0 

(V) 

CO  +  OH 

-^•C02  +H 

(8) 

Figure  12  Effects  of  y*0  on  the  temperature  rise  (ATmax  = 
Jmax  -  T*0 f )  over  the  whole  computational  domain 


Due  to  the  relatively  high  activation  energy  of  reaction 
(7)  (about  60.3  kJ-moL1),  OH  occurs  primarily  in  rela¬ 
tively  high  temperature  region.  In  the  conditions  of 
298.15  K  and  101.32  kPa,  its  reaction  heat  is  about 
+70.5  kJ-moL1  (endothermic  reaction)  and  the  reaction 
heat  of  reaction  (8)  is  about  -104.3  kJ-moL1  (exother¬ 
mic  reaction).  At  the  position  of  the  maximum  OH 
mass  fraction  in  the  mean  reaction  zone,  much  OH  is 
produced  by  reaction  (7).  Then,  OH  is  consumed  in 
the  downstream  and  releases  much  heat  due  to  the  fol¬ 
lowing  reaction  (8).  Thus,  rOH  always  peaks  slightly 

farther  than  rT  . 

1  max 


Figure  13  shows  the  radial  location  of  the  maxi¬ 
mum  temperature  (rT  )  and  OH  mass  fraction  ( r0H  ) 


3.4  Effects  of  the  hot  coflow  temperature  ( T*0 f  ) 
on  the  flame  structure 
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Figure  13  Radial  location  of  the  maximum  temperature 
and  the  maximum  OH  mass  fraction  for  Cases  1-7 
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Figures  14  and  15  show  that  as  7c*f  decreases, 
the  peaks  of  the  temperature,  OH  and  CO  mass  frac¬ 
tions  decrease  while  the  peak  of  the  H2CO  mass  frac¬ 
tion  increases.  That  is,  decreasing  Tc*f  generally  leads 
to  less  consumption  of  CH4,  thus  produces  less  OH 
and  CO.  Here  the  effect  of  7c*f  on  H2CO  is  contrary 
to  +o2 .  This  is  because  when  temperature  is  high, 
H2CO  is  consumed  rapidly  and  resides  quite  shorter, 
so  cannot  be  measured. 

Figures  16  and  17  show  that  Cases  8-12  should 
be  in  classical  mode.  For  Case  13  with  7c*f  =  1200  K, 
at  v=30  mm,  the  temperature  peak  is  only  slightly 
above  7c*f  (about  70  K)  and  a  very  small  amount  of 
OH  radicals  (about  2.5 xl0~4)  is  produced  due  to 
combustion;  but  at  v  =  200  mm,  the  temperature  rise  is 
about  400  K  and  OH  mass  fraction  is  about  6.0xl0"4. 
Thus,  Case  13  should  be  in  classical  mode. 

Figure  18  shows  the  effects  of  7c*f  on  the  tem¬ 
perature  rise  (  ATmax  =  Tmax  -  T*of )  over  the  whole 
computational  domain.  It  is  know  that  specific  heats  of 
flue  gases  such  as  C02,  N2,  02  and  H20  all  increase 
with  40f ,  so  that  the  same  amount  of  heat  releasing 
from  combustion  causes  the  temperature  to  rise  less 
when  the  coflow  temperature  is  higher.  At  the  same 
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Figure  14  Temperature  field,  OH  distribution  and  H2CO  distribution  in  the  whole  computational  domain  with  different  Jc*of 
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CH4  distribution  and  CO  distribution  in  the  whole  computational  domain  with  different 


time,  the  combustion  reaction  rate  increases  as  7c*f 
increases.  Whether  the  temperature  rise  (A7max)  over 
the  whole  domain  increases  or  decreases  depends  on 
which  effect  is  more  significant.  As  shown  in  Fig.  18, 
when  7c*f  is  higher,  A 7max  over  the  whole  domain 
tends  to  be  smaller;  in  other  words,  the  temperature 


distribution  tends  to  be  more  uniform.  Nevertheless, 
Fig.  18  indicates  that  in  industrial  applications,  the 
higher  the  coflow  air  is  preheated,  the  more  uniform 
the  temperature  distribution  is  achieved. 

Figure  19  shows  that  as  7c*f  decreases,  r0H 
and  rT  are  almost  to  be  the  same.  Varying  71*  f  has 
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Figure  16  Radial  distributions  of  the  mean  temperature 
(7)  for  Cases  8-13 
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Figure  17  Radial  distributions  of  the  mean  OH  concentra¬ 
tion  (voh)  for  Cases  8-13 
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little  effect  on  rnH 

UMmax 

also  shows  that  r0H 


and  rT  .  Similar  to  Fig.  13,  Fig.  19 

1  max 

is  slightly  larger  than  rT  . 

max  1  max 


4  CONCLUSIONS 


The  present  study  has  numerically  investigated 
the  effects  of  the  coflow  O2  level  ( )  and  temperature 
(7c*f)  on  non-premixed  MILD  combustion  produced 


by  a  JHC  burner  [15]  fed  with  the  CH4/H2  mixture. 
The  EDC  model  with  detailed  chemical  mechanism 
(modified  GRI-Mech  3.0)  is  implemented  for  the 
simulation.  The  computational  results  are  in  good 
agreement  with  the  experimental  data  [15].  Based  on 
the  analysis  on  the  computational  results,  several 
conclusions  can  be  drawn  below: 

(1)  As  y*0  decreases,  the  radial  profiles  of  tem¬ 
perature  distribution  is  more  uniform,  the  combustion 
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Figure  18  Effects  of  Jc*of  on  the  temperature  rise  ( A  Jmax  = 
Tmax  -  Jc*f )  over  the  whole  computational  domain 
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Figure  19  Radial  location  of  the  maximum  temperature 
and  the  OH  mass  fraction  for  Cases  8~13 
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mode  can  be  transformed  from  classical  combustion 
mode  to  MILD  mode;  the  peaks  of  temperature,  in¬ 
termediate  specie  OH,  H2CO  and  CO  mass  fraction  in 
the  computational  region  are  severely  suppressed.  It 
takes  a  longer  distance  to  bum  out  CH4  and  CO  when 
Jq  decreases.  If  jPq  is  high  enough  ( y^  >80%), 
increasing  y q  leads  to  little  combustion  temperature 
increase. 

(2)  As  7c*f  decreases,  the  peaks  of  temperature 
and  OH  mass  fraction  decrease  while  H2CO  mass 
fraction  peak  increases.  When  7c*f  is  higher,  the  tem¬ 
perature  distribution  over  the  whole  domain  trends  to 
be  more  uniform.  It  suggests  that  in  industrial  applica¬ 
tions,  the  higher  the  combustion  air  is  preheated,  the 
more  uniform  the  temperature  distribution  is  achieved. 


(3)  Reducing  y^  causes  the  radial  locations  of 
the  maximum  OH  mass  fraction  and  the  maximum 
temperature  further  away  from  the  axis,  while  varying 
7c*f  affects  them  quite  slightly.  It  is  interesting  that 
the  radial  location  of  the  maximum  OH  mass  fraction 
is  always  slightly  farther  away  from  the  axis  than  that 
of  the  maximum  temperature. 

NOMENCLATURE 


'OH 


rji* 

Tof 

T 

Tmax 

A  Tm 

X 

Toh 

•To, 


Tco 

* 

To2 

£ 


2  —2 

turbulence  kinetic  energy,  m  -s 

radial  location  of  the  maximum  OH  mass  fraction,  mm 
radial  location  of  the  maximum  temperature,  mm 

coflow  temperature,  K 
temperature,  K 

maximum  temperature  in  the  whole  computational  domain,  K 
temperature  rise  A Tmax  =  7max  -  T*of ,  K 
axial  distance  from  the  fuel  exit 
OH  mass  fraction 

02  mass  fraction 

CO  mass  fraction 

coflow  oxygen  mass  fraction 

dissipation  rate  of  turbulence  kinetic  energy,  m2-s~2 
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